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Abstract 

We develop a theory of resonant inelastic x-ray scattering (RIXS) at the K edge in La2Cu04 
on the basis of the Keldysh Green's function formalism. In our previous analysis (Phys. Rev. B 
71, 035110 (2005)), the scattering by the core-hole potential was treated within the Born approx- 
imation, and a crude-model density of states was used for the 4p band. We improve the analysis 
by taking account of the multiple scattering in Cu3d-02p bands and by using a realistic 4p DOS 
obtained from a band calculation. The multiple scattering effect is evaluated with the use of the 
time representation developed by Nozieres and De Dominicis. It is found that the multiple scat- 
tering effect makes the K-edge peak in the absorption coefficient shift to the lower energy region 
as a function of photon energy, that is, the photon energy required to excite the Is electron to 
the -ftT-edge peak reduces. It is also found that the multiple-scattering effect does not change the 
two-peak structure in the RIXS spectra but modifies slightly the shape as a function of energy 
loss. These findings suggests that the multiple scattering effect could mainly be included into a 
renormalization of the core-level energy and partly justify the Born approximation, leading to a 
future application to the RIXS in three-dimensional systems. 

PACS numbers: 78.70.Ck, 74.72.Dn, 78.20.Bh 
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I. INTRODUCTION 



Resonant inelastic x-ray scattering (RIXS), taking advantage of the high intensity of syn- 
chrotron sources, has become a powerful tool to probe charge excitations in solids.-'^i'^i^'^i^ 
The momentum dependence is detectable in transition-metal compounds by using the K- 
edge resonance, because the wavelength of photon is an order of crystal lattice spacing. The 
spectra consist of several peaks as a function of energy loss, and the peak positions some- 
times move with changing momenta transferred to the crystal. This is quite different from 
the optical conductivity,- where the momentum transfer is limited to nearly zero. Electron 
energy-loss spectroscopy can also detect the momentum dependence on charge excitations, 
but it suffers from strong multiple scattering effects.- Therefore, RIXS is quite valuable to 
investigate charge excitations, which is comparable to the neutron inelastic scattering for 
detecting spin excitations. 

In the RIXS process, the Is core electron is prompted to an empty 4p state by absorbing 
the photon, then charge excitations are created in order to screen the core-hole potential. 
Finally the photo-excited 4p electron is recombined with the core hole by emitting the 
photon. Charge excitations are left with energy and momentum transferred from the photon 
at the end. Most theoretical studies on the momentum dependence of the RIXS spectra 
in cuprates have been based on the numerical diagonalization method for small clusters. 
The 4p band was replaced by a single level there.— ^iii^^ A single-band Hubbard model 
has been sometimes used with replacing the charge-transfer band (so-called "Zhang-Rice" 
band) by the lower Hubbard band in two-dimensional cases.— It is not clear whether this 
replacement is appropriate, since the low energy peak is assigned to an excitation from the 
charge-transfer band to the upper Hubbard band. The single-band Hubbard model cannot 
describe another high-energy peak, either. 

By contrast, in our previous papers,— we have formulated the RIXS spectra on the 
multiband tight-binding model by adapting the resonant Raman theory developed by 
Nozieres and Abrahams. The formalism was based on the Hartree-Fock approximation 
(HFA) to describe the electronic states in the antiferromagnetic (AF) phase and the random- 
phase approximation (RPA) to the propagation of electron-hole pairs. The HFA is known 
to provide a good starting point for undoped materials. The calculation reinforced by the 
RPA correction has reproduced well the experimental RIXS spectra as a function of energy 
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loss and the dependence on momentum in La2Cu04.— The same formula has been applied 
to doped systems.— Thus, the formula seems promising to analyze the RIXS spectra in more 
complicated systems in three dimensions. But it should be noted that we have used a crude 
model for the Ap band and have treated the core-hole potential in the intermediate state by 
the Born approximation with neglecting the multiple scattering effect. Since the core-hole 
potential is not weak, higher-order effects might be important. Therefore, it is desirable 
to check and improve the approximations made before going to further applications. One 
purpose of this paper is to address this issue by using a realistic 4p density of states obtained 
from the band calculation and by evaluating the higher-order effects. Another purpose of 
this paper is to give a detailed derivation of the formula of RIXS spectra, because details 
have been sketchy in our previous papers.—*^ We present the formula with an emphasis of 
the time representation. 

We invent a numerical method to treat the multiple scattering by the core-hole potential 
along the line of Nozieres and De Dominicis,— deriving the Dyson equation to the one- 
electron Green's function with the core-hole potential working in a finite time interval. 
The point is that, different from metallic systems, we can numerically solve the Dyson 
equation for insulating systems, because the equation has no singular term. On the basis 
of the solution, we evaluate the multiple scattering effects on the RIXS spectra as well as 
the absorption coefficient. This type of analysis has not been attempted before. In the 
absorption coefficient, it is found that the i^'-edge peak moves to the lower energy region 
due to screening, and that the intensity at a higher energy region is a little enhanced as 
an antiscreening effect. When the incident photon energy is tuned at the i^'-edge peak, the 
RIXS spectra have peaks at around 2 and 5 eV as a function of energy loss. The 2eV- 
peak corresponds to the excitation of electron from the charge transfer band to the upper 
Hubbard band. It is found that the two-peak structure does not alter but the shape is 
slightly modified by the multiple scattering effect; the 2 eV peak is enhanced while the 5 eV 
peak is suppressed. The enhancement of the 2 eV peak is reasonable, because the core-hole 
potential is most likely to be screened by the excitation of electron from the charge-transfer 
band to the upper Hubbard band in the intermediate state when the incident photon energy 
is tuned at the K edge. However, the multiple scattering effect is limited within a small 
correction in the RIXS spectra. This may partly justify the Born approximation after a 
renormalization of the core-level energy. 
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The present paper is organized as follows. In Sec. [TTl we introduce a model and formulate 
the RIXS spectra. In Sec. Illlt the electronic structure is calculated on the d-p model within 
the HFA in the AF phase of La2Cu04. The RIXS spectra are calculated within the Born 
approximation in comparison with the experiment in Sec. IIV[ The multiple scattering effect 
is evaluated on the absorption and the RIXS spectra in Sec. El Section |Vl] is devoted to 
the concluding remarks. 



II. FORMULATION FOR RIXS SPECTRA 



A. Description of model 

We start by the expression of the Hamiltonian of photon, 

Hph = ^uj^cl^c^a, (2.1) 

where operator Cqa represents the annihilation operator of the photon with momentum 
q, polarization a. For the interaction between photon and matter, we consider the dipole 
transition at the K edge, where the Is core electron is excited to the Ap band with absorbing 
photon and the reverse process takes place. This process may be described by 

H^ = ^Y1 ^= E ei"Vl,.s,<,Cq,e^'i-'-^ + H.c, (2.2) 



qa V H jr]o- 

where e^T'' represents the 77th component (r/ = x, y, z) of two kinds of polarization vectors 
(a = 1,2) of photon. Since the Is state is so localized that the Is Ap dipole transition 
matrix element is well approximated as a constant w. Annihilation operators p'j^^ and Sja- 
are for states and state Is at Cu site j, respectively. 

The Hamiltonians for the core electron and the 4p electron are given by 

His = eu^s]^Sja, (2.3) 
H,, = 5^eI,(k)p'i^X..- (2.4) 

krjcr 

The photo-created Is core hole induces charge excitations through the attractive core-hole 
potential, which may be described by 

Hu-3d = VY1 4o^ios]^,Sj^/. (2.5) 

jacr' 



Here V may be 5 — 10 eV in La2Cu04. We neglect the interaction between the core-hole and 
the Ap electron, since the 4p states are well extended in space with the bandwidth as large 
as ~ 20 eV. The excited 4p electron is finally recombined with the core-hole with emitting 
photon. In the end, charge excitations remain with receiving a momentum and an energy 
from scattering photons. 

In La2Cu04, Cu and O atoms form a two-dimensional network shown in Fig. [H There 
is nominally one hole per Cu atom. To describe this situation, we consider only the x"^ — y"^ 
orbital at the Cu site, which hybridizes the cr-bonding 2p orbitals at O sites {^^d-p" model). 
The corresponding Hamiltonian may be expressed as 

Hdp = ^d^d]jj, + ep^p\^p^„ 

+ ^ tjid^j^pe^ + }l.c. + ^ tu>pl^pii„ 

+ UdY,d]^di^d]idn^ (2-6) 

where dj^, is the annihilation operator for the —y^ orbital with spin a at Cu site j, and p^^ 
is the annihilation operator of the 1p orbital with spin o at site i. For La2Cu04, the transfer 
energies between and 2p orbitals and that between Ip orbitals are estimated from 
the local density approximation (LDA) calculations, that is, td = 1-3 eV, tp = 0.65 eV.-- We 
assume Ud = 11 eV, and the 3d level energy relative to the 2p level -|- \Udnd — = —0.7 
eV with Ud = J2a(d]adj<^) ■ These values are set the same as those in our previous paper.— 



B. Keldysh-Schwinger formalism 

Following Nozieres and Abraham,- - we use the Keldysh-Schwinger formalism to the RIXS 
spectra. We prepare the initial state that one photon exists with q^, in addition to a 
material in the ground state, which may be expressed as 

m=clJO), (2.7) 

where |0) is the ground state of the matter with no photon. Let H = Hph + Hdp + His + H^p 
be the unperturbed Hamiltonian of the system and H^ be the perturbation. Then the 5* 
matrix is given by 

U{t, -oo) =Texp!^-i j H^{t')dt'^, (2.8) 
5 




FIG. 1: Schematic view of the unit cell in the antiferromagnetic phase. Six orbitals (numbered 
1-6) are considered. There are two nonequivalent sites of Cu shown as A and B. The gray parts in 
orbitals indicate that their wave functions take negative values. 

with Hx{t) = exp(iHt)Hrcexp{—iHt). The probability of finding a photon with momentum 
q/, polarization at time to is given by 

^q,a,;q,a.(to) = f/(-00, to)<^,„, Cq,a, f/(to , "Oo) | $i) " (2-9) 

Expanding the S-matrix to second order in H^, 

U{t, -oo) = l + I H^{t')dt' 

J — oo 

+ f f T{H,{t')H,{t"))dt'dt", (2.10) 

we insert this into Eq. (12.91) . Figure [2] shows a schematic representation of the expansion, 
where the wavy lines indicate photons, which carry energy and momentum, and the solid 
lines with "4p" and "Is" represent the Green's functions of the 4p electron and the Is core 
hole, respectively. The upper and lower halves of the graph correspond to the so called 
"outward" and "backward" time legs, respectively. By factoring out the dependence on the 
photon frequencies, we obtain 

/to pu pto pu' 

du dt du' dt'5(t,M;tV)e^->(*'-*)e-^-/("'-"). (2.11) 
oo J —oo J —oo J —oo 
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FIG. 2: Expansion of the S matrix. The wavy hnes represent photon Green's functions. The sohd 
hnes with 4p and Is represent Green's functions of the 4p electron and the Is core hole, respectively. 

The transition probability per unit time with ^ oo is given by fixing one time, for 
instance, u = 0, 

/O roc pu' 

dt du' dt'5(t,0;tV)e*^'(*'-*)e-*"/"'. (2.12) 
oo J —oo J —oo 

III. HARTREE-FOCK APPROXIMATION 

The undoped cuprates show the AF long-range order at T = 0. It is known that the HFA 
works rather well in the AF phase even for the case of large Coulomb interaction between 
3(i electrons. Since the following analysis of RIXS is based on the HFA, we summarize the 
electronic structure within the HFA in this section. The unit cell in the AF phase contains 
six orbitals for each spin, as shown in Fig. [H Introducing the Fourier transform in the 
magnetic Brillouin zone (BZ), we define the single-particle Green's function for six orbitals 
in a matrix form, 

[G"^(k, u)],, = I (T(Ak,.(t)AL(0)))e-Mt, (3.1) 

where Ak^io- represents the annihilation operator of the electron with the orbital at site /i in 
the unit cell. For example, at Cu sites (/i = 1 and 4), it is given by 
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with i running over the sites "/i" of N/2 unit cells. 

Applying the equation-of-motion method to the Green's functions we obtain the relation, 

{ooi - J"(k))G"(k, io) = J, (3.3) 

where / is the unit matrix, and J'^(k) is given by 

ep -2tpCos(^-^) -tde''^ 2tpCos(^ + ^) 

ep tde''-^ 2t„cos(% + ^) 



ed + UdU^^ tdd'"^ -tde' 



2 



ep -2tpCos(^-^; 



(3.4) 



with 



1^ / |G"(k..)]„e-*^. = 1^ / P(k..)]..e-*^^. (3.5) 

Since = n^^, we may put 



k 



Lower triangle components are the Hermitian conjugates to the upper triangle components, 
which are omitted in Eq. (13. 4p . The J'^(k) is diagonalized by an unitary matrix ?7°"(k), that 
is, [U~^JU]jjf = EjllijSjji. Then the Green's function is expressed as 

G"(k, u) = t/"(k)D(k, u})U''{k)-\ (3.7) 

with 

Figure [3] shows the dispersion relation -E'j(k) as a function of k along symmetry lines. The 
conduction band and the top of the valence band may be called the "upper Hubbard" band 
and the "charge-transfer" band, respectively. 



IV. BORN APPROXIMATION FOR RIXS SPECTRA 



We consider the process that an electron-hole pair is singly created to screen the core-hole 
potential. The multiple-scattering effect beyond the Born approximation is neglected. Figure 
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FIG. 3: Dispersion relation of Cu3d-02p bands along symmetry lines in the Brillouin zone. 

m displays the corresponding diagrams, where the shaded part represents the Keldysh-type 
Green's function, 



-iuj(s'-s) 



(Pq/.'a'(s')P-qM<x(s)), 



27r 



with 



k+q/x(T"'k/X( 



(4.1) 
(4.2) 

(4.3) 



The momentum conservation requires the relation q = q?; — q/, and k runs over the mag- 
netic first BZ. The superscripts + and — stand for the backward and outward time legs, 
respectively.— In the lowest order, it is given by 



k 

and thereby 

>?;,^l(q,^) = 5..'$^$^5(^-i?,([k + q]) + i?,(k))[l-n,,([k + q])K-(k) 

X u;,^, ( [k + q] ) ij;;, ( [k + q] ) u;^ (k) f/;* (k) , 



(4.4) 



(4.5) 



where j and f stand for energy eigenstates, and [k + q] is the reduced value of k + q into 
the magnetic first BZ by a reciprocal lattice vector G, that is, k + q = [k + q] + G. Since 
two Cu sites in the unit cell are involved to excite an electron-hole pair, phase factor e**^** is 
attached to the wave function at the B site. 



[/;^.([k + q])=e^^'^f/;,.([k + q]) for /x = 4, 
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(4.6) 




s 




f 4p u 



FIG. 4: Diagrams for the RIXS intensity within the Born approximation. The dotted lines represent 
the core-hole potential. The shaded part represents the Keldysh-type Green's function, which 
connects the outward time leg on the top half and the backward time leg on the bottom half. 

where /j, — A corresponds to the Cu B site, and a [= (a, 0)] represents a position of the B 
site relative to the A site in the unit cell. 

The product of Green's functions of the Ap electron and the core hole gives simply a factor 
exp[i(e4p(p) — eis — iFig — coijt] on the outward time leg and a factor exp[— i(e4p(p) — ei^ + 
iVis — uJi}{t' — u')] on the backward time leg, where Fi^ is a lifetime broadening width of 
the Is core hole. In a conventional calculation based on the Fourier transforms, this feature 
is hidden. This property will be fully utilized in the study of the multiple scattering effect 
in the next section. The Keldysh Green's function carries the time dependent factor e*'^* to 
the outward time leg and e^*'^*' to the backward time leg. Note that the core-hole potential 
works only in intervals [t, 0] and \t\u']. Integrating the time factor combined to the above 
product of Green's functions, with respect to s and t in the region of i < s < 0, — oo < t < 0, 
we obtain 



Here the sum over Ap states is replaced by the integration of the 4p DOS projected onto the 




(4.7) 
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t] (= X, y, z) symmetry, plp{e). A similar factor has been derived in third-order perturbation 
theory by Abbamonte et al.^° The integration with respect to s' and t' in the backward 
time leg gives the term complex conjugate to Eq. fl4.7p . The integration with respect to u' 
gives the energy-conservation factor, which guarantees that uj in Eq. (14 .yp is the energy loss, 

00 = UJi — iOf. 

We evaluate the higher-order effect on ^^t^^o-(Q') = (^)^)] within the RPA. Figure [5] 
shows the corresponding diagram. Collecting up the ladder diagrams, we obtain 



where 



with 



A 



Aiia-1,M20"2 * 



(g)= I-UdF{q) 



Mia"l,M2a"2 



(Fl{q) FlM \ 
FUq) 

FlM Fl{q) 

V F^q) FUq)) 

f Ud \ 
f/d 



f/d 

V f/d / 

Here the polarization propagator F^^, (g) is given by 



N E •(k)^,V([k + q])f/,7,([k + q]) 



n,(k)[l-n,,([k + q])] 



n,,([k + q])[l-n,(k)] 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



u - Ef{[k + q]) + Ej{k) +i6 uj - Ey{\k + q]) + ^^(k) - id _ 
Combining all these relations together, we finally obtain an expression of the RIXS in- 
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FIG. 5: Vertex function within the RPA; (a) the channel of keeping spin, and (b) the channel 
of reversing spin. Broken lines between bubbles represent the 3d Coulomb interaction, which is 
effective only between 3d electrons with opposite spins. 

tensity, 

W{qi,a,;qf,af) = 27r V V V + i?,(k) - i?,,([k + q]))n,(k)[l - n,,([k + q])] 

X A;^^,^,,,(g)f/;;.(k)f/;;^.(k)f/;^^.([k + q])f/;^([k + q])A^,^,^.(g) 



X 



(4.13) 



where the incident photon has the momentum and energy = (qj,C(jj), polarization e'-"'^ 
and the scattered photon has the momentum and energy qj = (qj,u;/), polarization e^°'f\ 
The transfer of momentum and energy is denoted as q = q^ — q^, uj = uji — uj f . 
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In 



We have already reported the same formula without a detailed derivation in Ref. 
that study, we calculated the RIXS spectra of La2Cu04 by assuming the two-dimensional 
cosine dispersion for the Ap band and tuning the incident photon energy to excite the core 
electron to the peak of the Ap DOS. We improve the treatment by using the realistic 4p DOS 
calculated with the LDA. Although the LDA fails to predict the AF insulating ground state, 
the 4p band is expected to be properly described because of weak correlation in wide bands. 
Figure [6] displays the 4p DOS convoluted by the Lorentzian function with the half maximum 
full width (HMFW) 2ris = 1.6 eV. Reflecting the tetragonal structure, the p^-symmetric 
DOS is found quite different from the p^- and Pj^-symmetric DOS's. Since the experiment 
by Kim et al. has been performed with the polarization along the z axis, we concentrate our 
calculation to the p^-symmetric DOS. Thin lines represent the absorption coefficient with 
taking account of the Is — » 4p dipole matrix element calculated by the band calculation 
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Energy (eV) 



FIG. 6: Density of states for the 4p band convoluted with the Lorentzian function with HMFW 
2ris = 1.6 eV. The origin of energy corresponds to the bottom of the 4p band. Thin hues represent 
the absorption coefficient with taking account of the dipole matrix element evaluated by the band 
calculation. 

(the origin of energy corresponds to the excitation of the Is electron to the bottom of the Ap 
band). The curves are very close to the 4p DOS, indicating that replacing the dipole matrix 
element by a constant value in Eq. (12.21) is a good approximation. 

We calculate the RIXS intensity from Eq. (I4.13P by tuning the incident photon energy 
at the peak of the pz DOS, the 6 function is replaced by the Lorentzian function with 
HMFW= 0.2 eV in order to take account of the instrumental resolution. Figure [7| shows 
the calculated spectra as a function of energy loss for momentum transfer (0, 0), (tt, 0), and 
(vr, vr). It should be noted here that changing magnitude of V does not change the spectral 
shape within the Born approximation; as is evident from Eq. (I4.13p . W{qi,ai;qf,af)/V'^ is 
independent of V. The present result using the realistic Ap DOS is nearly the same as our 
previous one^ using a two-dimensional cosine dispersion for the 4p band. We obtain the 
spectra consisting of two peaks at about 2 and 5 eV as a function of energy loss, whose shapes 
move with changing momenta in good agreement with the experiment. The 2 eV peak is 
assigned to the excitation from the charge-transfer band to the upper Hubbard band. The 
2eV peak at (0, 0) is considerably enhanced, while that at (vr, vr) is considerably suppressed 
by the vertex function of RPA in agreement with the experiment. As already emphasized 
in our previous paper,— the intersite correlation reinforced by the RPA is quite important 
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5 10 

Energy Loss (eV) 

FIG. 7: RIXS spectra as a function of energy loss uj = LOi — i^Jf, within the Born approximation. 
Momentum transfer is (0,0), (tt, 0), and (tt, tt). The incident photon has the z-polarization and 
the energy corresponding to the excitation from the Is core to the peak of the Pz DOS. The dotted 
and sohd hues are the results of the HFA and RPA, respectively. 

for the quantitative understanding of the RIXS spectra. The difference between the present 
and the previous results is that the 2 eV peak is a little enhanced relative to the 5 eV peak 
in the present result. 
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V. MULTIPLE SCATTERING BY THE CORE-HOLE POTENTIAL 



Since the Is core-hole potential is rather strong, many electron-hole pairs may be created 
in the Cu3c?-02p bands in order to screen the potential. This process may influence not only 
the absorption spectra but also the RIXS spectra. 

It is known in metallic systems that the creation of an infinite number of electron-hole 
pairs leads to "singular" behavior in the absorption spectra.— In insulating systems, how- 
ever, the number of electron-hole pairs is limited due to the finite energy gap, and the sin- 
gularity would not come out. Under this condition, we calculate the RIXS and absorption 
spectra on the basis of the time representation developed by Nozieres and De Dominicis.— 

Since the potential is localized at Cu sites, it is convenient to introduce the one-electron 
Green's function for the 3d electron localized at the /ith site (yU = 1 and 4), 

f^,a{s2,Sl) = -i{T{df,„{s2)dl^{si)). (5.1) 

Here the core-hole potential works only between the time interval [t,0], and t < si, S2 < 0. 
Figure [H] shows a diagram contributing to '^^a{.S2, Si). Collecting similar diagrams, we obtain 
the Dyson equation 

^M<x(52,si) = G(fj(s2-si)-\/ G^;}{s2-S3)ip^^{s3,s,)dss, (5.2) 
where G^^}{s) represents the unperturbed local Green's function, which is given by 

G^^Hs) = (-^)|Elf^.^•(k)l'(l-^.(k))e-^^^^'^^ ^>0, (5.3) 

ki 

= 4 5^ |/7;,■(k)|S(k)e-^^^«^ s < 0. (5.4) 

kj 

In metallic systems, G^/^}{s) behaves as 1/s for |s| >> D^^ with D being the bandwidth. 
In this situation, the Dyson equation becomes a singular integral equation, which needs 
special care.— On the other hand, in an insulating system we are dealing with, the asymptotic 
term of 1/s would vanish for s >> d^^ > D^^ with d being the band gap, owing to the 
presence of the energy gap. In addition, the solutions for large values of |t| are less important 
in the absorption coefficient and the RIXS spectra, since the finite core-hole life time reduces 
their contributions. Therefore, we are allowed to solve the integral equation [Eq. (15.21) ] 
in the finite interval [t, 0] with setting the maximum value of |t| by 12 (eV)~^ In the 
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FIG. 8: Diagram for the one-electron Green's function 99'^(s2,si). Solid circles represent the core 
hole potential working on the electron. 

numerical solution, we divide the interval into small widths by 0.02 (eV)^^. The reduction 
factor at \t\ = 12 (eV)~^ due to the core-hole life-time width Fi^ = 0.8 eV is estimated as 
exp(— 12Fis) = 0.000067, which seems sufficiently small. Note that we need the one-electron 
Green's functions only for S2 = t or Si = t in later use. 

The appropriate value of V is not known and may be model dependent. It may be better 
to take the value as a parameter. In the following, we calculate the absorption and RIXS 
spectra for V = 6 and 9 eV, in order to get some insight into the parameter dependence. 

A. Absorption coefficient 

Before going to the RIXS spectra, we first formulate the absorption spectra. Let a photon 
have frequency Ui and polarization direction rj. Then, the absorption coefficient is calculated 
from 



Figure shows a corresponding diagram, where the interaction between the core hole and the 
4p electron is neglected. To screen the core-hole potential, electron-hole pairs are created. 




(5.5) 
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FIG. 9: Diagram for the absorption coefficient at the Cu K-edge. A photon enters with frequency 
CO and polarization rj. The soHd fines with 4p and Is represent the Green's function of the 4p 
electron and Is core hole, respectively. Other solid lines represent tfie local Green's function. The 
dotted lines denote the core-hole potential. 

Without such bubble diagrams, the absorption coefficient is simply given by 
A^{uJi) = -- — Re 22 / exp{i(e2p(k) - ei^ - iTu - uji)t}dt, 

-de. (5.6) 



\w?Vu 



P4p(e) 



27ru;j J {uji-t + eu)"^ + Tl^ 

This is nothing but the 4p DOS convoluted with the Lorentzian function with HMFW ITis. 

The electron-hole pairs are treated on the basis of the linked cluster expansion. The 
contribution is simply given by a factor e*" , where 

C{t) = -\v'y. f [\,As2,si)G^Si',-S2)ds,ds2. (5.7) 



In the actual calculation of C{t), we numerically evaluate dC{t)/dt from the relation 
dC{t) 



dt 



E / [V,.{t, s)G'^^{s -t) + ^,^{s, t)G<-^\t - s)] ds, (5.8) 



and then integrate this quantity with respect to t, instead of directly evaluating Eq. (15.71) . 
The product of Green's functions of the Ap electron and the core hole gives a factor 
exp[i(e4p(k) — eis — iTu — uji)t], which is independent of times at which dotted lines are 
attached to the core-hole line in Fig. O Therefore, the absorption coefficient is given by 

|2 

Urn Re[K''^{ujf,T)], (5.9) 



Ar,{uJi) 



\W\ 
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with 

K\uJuT) = J2 rexp[2(e2p(k)-ei,-zri,-a;,)t]e^Wdt 
k ^'^ 

= j p%{e)dej^ exp[^(e - ei, - iT^ - ^0^]e^Wdt. (5.10) 

We integrate numerically Eq. f lS.lOp up to |T| = 12{eV)^^ . As already mentioned, the 
reduction factor due to the core-hole lifetime is estimated as exp(— ris|T|) = 0.000067. Fig- 
ure (TU] shows the calculated absorption coefficient. The broken lines represent the quantity 
without taking account of the screening by electron-hole pairs [C{t) = 1 in Eq. flS.lOp . 
The curves reproduce well the curves with |T| ^ oo shown in Fig. [6l indicating that the 
introduction of the cutoff |T| = 12{eV)~^ causes only minor errors. The thick and thin 
solid lines represent the absorption coefficient for V = Q eV and 9 eV. The K-edge peak 
moves about 1.5 eV for V = Q eV and about 3 eV for V = 9 eV to lower-energy positions, 
which may be called "well- screened" peaks. In addition, other peaks appear at about an 8 
eV-higher position for V = Q eV and at about an 11 eV-higher position for V = 9 eV, in the 
2;-polarization curves, probably owing to antiscreening. These peaks may be called "poorly 
screened" ones. However, the poorly screened peaks are weak with broad widths and nearly 
merged into the background. Therefore, the screening effect may be, on the whole, taken 
into account by a renormalization of the core-level energy. 



B. RIXS intensity 

Confining ourselves within the HFA, we study the multiple scattering effect on the RIXS 
spectra in this subsection. Figure [11] shows the effect of the core-hole potential. Both the 
Born and the multiple scattering terms accompany the excitation of electron-hole pairs. The 
lines with Ei and E2 denote the factors e~*-^^*^ and e'^^'*^ coming from the Keldysh-type 
Green's function, ^^Xr'*'°l(5')5 which combines the outward and backward time legs. The 
product of Green's functions of the 4p electron and the core hole, which is hidden in the 
background, gives a factor exp[z(e4p(k) — eis — iVig — uji)t\. The electron-hole pairs give 
an extra factor e*"^*-*, which is the same as in the absorption coefficient. Combining these 
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FIG. 10: Absorption coefficient for z, and x, y polarizations. The broken line represents the result 
without the screening. The thick and thin solid lines represent the results with the screening by 
electron-hole pairs for V = 6 and 9 eV, respectively. The origin of energy corresponds to the bare 
energy exciting the Is electron to the bottom of the Ap band, that is, e^p — eis with e^p being the 
energy of the 4p band bottom. 

factors, we obtain 



Jt Jt Jt 

[K^{uJi;-oo)- K^{uJi;t)]dt 

o 



S2 I e'^^^V"(s2,si)e- 



(5.11) 



X |i/e*(£^2-^i)* -V^ [e'^^'Vls, t)e-''''' + e'^^ V"(i, s)e-'^i^] dsj , (5.12) 

where K^{uji;t) is defined by Eq. (15.101) . The multiple scattering effect is included in 
V'^{s2,Si). We finally obtain the formula of the RIXS intensity with replacing L^^{uJi;uj) 
by L^{uf,E,{k),Ey{[k + q])) in Eq. (glS]). 

On the basis of this formula, we calculate the RIXS spectra. We set the core-hole potential 
to be y = 6 and 9 eV, and tune the incident photon energies to correspond to the peak in 
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(a) 

t 




FIG. 11: Multiple scattering process for RIXS;(a)tlie Born term accompanying electron-hole pairs, 
and (b)a multiple scattering term. The solid lines represent the one-electron Green's function Gq 
and solid circles represent the core-hole potential causing scattering. 

the absorption coefficient shown in Fig. [T0|, which are shifted from the unscreened positions. 
Figure [12] shows the calculated spectra divided by V^. Dividing by V"^ makes the spectra of 
the Born approximation independent of the value of V . 

The 2 eV peak is enhanced and the 5 eV peak is suppressed by the multiple scattering. 
This tendency becomes stronger with increasing values of V . This may be related to the 
fact that the incident photon energy is tuned to the well-screened peak, where the excitation 
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from the charge-transfer band to the upper Hubbard band is most effective for screening. 
Nevertheless, the correction remains minor on the RIXS spectral shape. By comparison with 
the result in the previous section, the RPA correction is found more important. 

VI. CONCLUDING REMARKS 

We have developed a theory of the RIXS on the basis of the Keldysh Green's function 
formalism. In this formula, we have described the electronic structure by the d-p model 
within the HFA, and have treated the electron correlation within the RPA. The RIXS 
spectra can be interpreted as a band-to-band transition reinforced by the RPA. We have 
analyzed the RIXS spectra of La2Cu04 using the realistic 4p DOS, in good agreement with 
the experiment. 

One of the important approximations in our theory was the use of the Born approximation 
to treat the core-hole potential in the intermediate state. The multiple scattering effect may 
not be small, since the core-hole potential is by no means small. For evaluating the multiple- 
scattering effect, we have invented a numerical method along the line of Nozieres and De 
Dominicis. In the absorption coefficient, we have found that the fT-edge peak moves to a 
lower-energy region due to screening (well-screened peak), and that another peak appears 
at a higher-energy region due to antiscreening (poorly screened peak). However, the latter 
peak is weak with a broad width and nearly merged into the background. On the RIXS 
spectra, we have found that the multiple scattering modifies the spectral shape slightly. 
These findings suggest that the multiple scattering effect could be mainly included into a 
renormalization of the core- level energy, and partly justify the use of the Born approximation. 
This paper is the first attempt to evaluate multiple scattering effects. The above conclusion 
is important, since the Born approximation is employed without examining the validity in 
several attempts of the RIXS analysis. The Born approximation is expected to work well 
in other insulating transition- metal compounds, such as NiO, LaMnOs, etc. The present 
formula seems promising to analyze the RIXS spectra for such systems, because it can deal 
with multiband tight-binding models in three dimensions. 

The HFA and RPA are known to work well in the presence of the AF long-range order. 
Upon doping, the AF order is easily destroyed in cuprates. Recently the doping effects on 
the RIXS spectra have been observed in cuprates.— "^"^ In such a situation, the HFA no 
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FIG. 12: RIXS spectra as a function of energy loss for momentum transfer (0, 0), (vr, 0), and (vr, vr), 
with the inclusion of the multiple scattering effect. The RPA correction is not included. The 
thick and thin lines represent the intensities divided by V'^ for V = 6 and 9 eV, respectively. The 
incident photon is assumed to have the z polarization and the energy corresponding to the peak 
in the absorption spectra shown in Fig. [TOj The broken lines represent the intensities divided by 
V'^ within the Born approximation, where the incident energy is set to be the unscreened K-edge 
peak. Dividing the intensities by V'^ makes the results of the Born approximation independent of 
the value of V. 
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longer provides a good starting point, and we have to go beyond the HFA-RPA scheme. A 
study along this line is now under progress. 
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